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t-h We introduce a family of glassy models having a parameter, playing the role of an 

interaction range, that may be varied continuously to go from a system of particles in 

CN 

d dimensions to a mean-field version of it. The mean- field limit is exactly described 

Qh 

<^ by equations conceptually close, but different from, the Mode-Coupling equations. 

We obtain these by a dynamic virial construction. Quite surprisingly we observe that 
in three dimensions, the mean-field behavior is closely followed for ranges as small 
as one interparticle distance, and still qualitatively for smaller distances. For the 

<D 

original particle model, we expect the present mean-field theory to become, unlike the 
Mode-Coupling equations, an increasingly good approximation at higher dimensions. 
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In the past few years, there has been considerable activity on the application of Mode- 
Coupling theory to liquid systems. In its original conception, Mode Coupling is an ap- 
proximation for the dynamics in which an (infinite) subset of corrections coming from non- 
linearities is taken into account. The theory has become popular not so much for the accuracy 
of its predictions - numerical confirmation often demands considerably good will to accept 
- but because it gives a unified and qualitative view of the first steps of the slowing down 
of dynamics, as the system approaches the glass transition. 

Mode-Coupling theory was originally confined to the dynamics in equilibrium liquid 
phase. However, similar approximations may be applied to the equilibrium statistical me- 
chanics of systems above and below the glass transition, and to the non-equilibrium ("ag- 
ing") dynamics below the glass transition. Elaborating on an idea of Kraichnarr , Kirkpatick, 
Thirumalai and Wolyned^ noted that one may view these approximate theories as being 
the exact description of the properties of an auxiliary model, different from the original one. 
These turn out to be disordered models having nature that is 'mean-field' in the following 
sense: if a particle (or spin) A interacts strongly with two other particles B and C, then 
B and C do not interact strongly with one another. This may happen either because the 
network of interactions is tree-like, or because all individual interactions are weak. 

There is a large family of such 'mean-field' models, going beyond the one that leads to the 
original mode-coupling equations, and they may all be treated with the tools developed in 
the context of spin glass theory. When applied to the glass transition, the whole strategy is 
referred to as "Random First Order" scenario: a single name for approximation schemes that 
may be different is justified because the expectation is that the nature of the glass transition, 
of the equilibrium glass phase, and of the out of equilibrium dynamics, is qualitatively the 
same in all these mean-field models. The wider question whether this scenario holds strictly 
at finite dimensions is still very far from established. 

Once one recognizes that Mode-Coupling is a form of mean-field theory, the first instinct 
is to ask under which conditions it becomes exact, in particular if it does so in high dimen- 
sions. The answer for the latter question is that it does notP. This lack of control over the 
approximation is problematic, because there is no unambiguous way of relating features of 
a realistic system with those of the Mode-Coupling solution - and we often are not sure 
whether some qualitative crossover in the behavior of an experimental system should be 
associated with the idealized Mode-Coupling transition, and in what sense. 
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In this paper we study an approximation of the same general mean-field class than, but 
different from, the one leading to the mode-coupling equations. In order to bridge the gap 
between this limit and reality, we build explicit models where an interaction range is tuned by 
some parameter, thus allowing to go continuously from mean-field to true finite-dimensions 
by varying this parameter. Work in this direction already exist for spin glasses where 
one can consider models with interactions with tunable range, as originally proposed by 
Kac. As we shall mention below, for particle systems, the usual program a la Kac meets a 
problem as the interactions are made longer in range and less strong: at low temperatures 
and large densities particles tend to arrange themselves in clusters^^, themselves arranged 
in a crystalline or amorphous 'mesophase' structure. Thus, it seems that in order to prevent 
this, one is forced to add a short-range hard-core repulsion, thus spoiling the Kac (mean- 
field) nature of the model. In this paper we follow a different path, based on a suggestion 
already made by Kraichnan fifty years ago: we study particles with short-range interactions 
which are, however, 'shifted' by a random amount having a typical range, which is our 
parameter. 

Within this framework we are able to address several issues related to the glass transition, 
revisiting them via the mean-field model we introduce. As an example, we are able to answer 
questions such as: "what is the relation between the point at which the dynamics becomes 
nonexponential and the dynamic transition", because we can continuously take the model 
from finite dimensional to a mean-field limit, a situation where both transition points are 
well-defined, independently of any fitting procedure. 

This paper is divided in two parts, analytic and numeric, which may be read indepen- 
dently. Sections II and III are devoted to the analytic treatment of the statics and the 
dynamics of the liquid phase, respectively. The main new result is an equation for the dy- 
namics that plays the role of the mode-coupling equation, and is exact in the mean-field 
limit. Sections IV and V present the numerical tests for statics and dynamics, respectively. 
We are able to compare the results in the mean-field limit with the ones for finite parameter 
A, all the way down to the ordinary particle model A = 0. In section VI we discuss an instance 
where having an approximation with some limit in which it is well controlled is reassuring: 
there has been some doubt whether the so-called "onset temperature" (or pressure at 
which the equilibrium dynamics becomes nonexponential (and the inherent structures start 
to have deep energies), should be identified with the mode-coupling transition. Here, by 
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taking continuously the parameter A to infinity, we find that they are in fact two distinct 
pressures - at least in the limit in which they are both well defined. 

I. MODEL 

A. Kac models, clustering and the Kirkwood instability 

One can introduce a mean-field treatment of particle systems in different ways. One may, 
for example, use explicit infinite range interactions^^. Glassiness is obtained by choosing 
a potential imposing a strong frustration, and there is in principle no need for quenched 
disorder. Next, one may consider long, but finite ranges, in the spirit of Kac interactions. 
Although quite intuitive, the choice of the potential is in practice difficult, as crystallization^ 
or instabilities in the liquid phase (like the Kirkwood instability ^ 12 * 17 ^ ) easily set in as 
soon as the interaction range is finite. 

Consider, for example, the model studied by Dotsenkd 14 ' 15 ! Particles are in a confining 
potential V^ on /(x a ) (which may be harmonic) and interact with a long-range, oscillatory 
potential: 

H = ^=X! COS (l X a- X b|) + E^on/(x a ) (1) 

V JS a*b a 

Dotsenko showed that the system indeed has a mean-field glass phase, induced by the frus- 
tration due to the conflict of attraction and repulsion. The next step, in an ordinary Kac 
program, would be to introduce the model with a finite range 7: 

H = -7= E COS d X a " X b |)e-I x - X ^ + Y, Kon/(x a ) (2) 

The result is disappointing: instead of giving a glass, the system arranges as follows: it 
forms clusters of many particles, themselves disposed in a crystalline arrangement, at the 
optimal distance so that the interaction between clusters is minimized. The way to avoid the 
crystal-of-clusters mesophase is, of course, to add a hard-core that hampers the clustering, 
but then the model is no longer mean field. Another related difficulty with this strategy is 
that even the liquid phase may have a transition to a 'liquid' with spatial modulation, the 
Kirkwood instability ^ 12 ' 17 ! This phase is not without interest of its own (see Fig. [3] for a 
numerical simulation), but it is not what we are wishing to study here. 

Another way to construct a mean-field model is to work with particles on a Bethe 
lattice®^. This introduces quenched disorder and the mean-field nature at the same time. 
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By increasing the graph connectivity, one can, at least formally, recover the original finite 
dimensional model by setting the graph connectivity to infinity. 

B. Kraichnan's proposal and beyond. 

In this paper, we will follow another route. We study family of models which are defined 
through the Hamiltonian: 

#({x},{A})= £ ^(xi-xj-Ay) (3) 

<i,3> 

where V is a short-ranged interaction potential. The x's are the positions of the particles 
and the A's are quenched random variables with a probability distribution P(|A|), which 
has a variance A 2 . We also impose that A^- = Ajj. The model for A = oo was first introduced 
by Kraichnan some fifty years ag (P. 

When A = 0, P(|A|) = <5([A|), and the model reduces to an usual ci-dimensional system. 
On the other hand, when A -> oo, one particle i can interact with particles j which are 
possibly anywhere in the system, as long as |xi - Xj - A^l is of the order of the range of 
the potential V. Therefore, the system tends to have a mean-field nature in this limit even 
though one particle effectively interacts with a finite number of other particles, as it does 
in a conventional finite-dimensional hard sphere system (see Fig. [T|. Thus by tuning A, the 
model Q goes from a finite dimensional system (A = 0) to a mean-field realization of the 
same system (A -> oo). 

In the liquid phase, the mean-field limit of the model has an entropy of an ideal gas plus 
only the first virial correction, just like a van der Waals gas. Physically, this comes from the 
fact that it is very unlikely that three (or more) spheres effectively interact simultaneously 
with one another, as is sketched in Fig. [2] For instance, in the hard-sphere case, one would 
need to have, for particles i, j, and k (having a diameter D): 



|Xi 






- D 


l x r 


-x k 


-Ajkl 


~ D 


|x k 


"Xi 


-Awl 





which requires that the random shifts satisfy |Ay + Aj k + A k j| ~ D. This of course is very 
unlikely for shifts |A| ~ L, where L is the linear size of the system. 
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FIG. 1. The black particle effectively interacts with a finite number of particles (in red) that have 
appropriate random shifts (blue arrows), but these particles may be anywhere in the sample. 




FIG. 2. Left : Effective three-body correlation in a finite-dimensional hard-sphere system (with 
no random shifts), leading to corrections beyond the first virial term. Center : Same situation 
in the random shift model, from the point of view of particle 1. Particle 1 sees particles 2 and 3 
nearby. Right : From the point of view of particle 2, there is no effective three-body interaction, 
as the random shifts displaces particle 3 very far from particle 2. 

In this work, we will concentrate mostly on the hard-sphere potential, but the procedure 
is a very generic way to obtain a mean-field limit, and can even be generalized to objects with 
rotational degrees of freedom, where one can introduce a rotational disorder. We studied 
both a monodisperse system, for conceptual simplicity; and a bidisperse one, to be able to 
work with arbitrary small A without having to deal with crystallisation. 

From now on, we set the diameter D of a sphere (in the monodisperse case) to D = 1, 



which will be used as length scale. The temperature T is set to 1 almost everywhere, as 
it is irrelevant for hard spheres. The only part where we keep it explicit is in the section 
dealing with the dynamics of the model. For simplicity we introduce the following notations: 
a sphere with diameter 1 in dimension d has a volume 2~ d Vd, and a surface sa- 
il. STATICS OF THE LIQUID PHASE 
A. Grand-canonical formalism and Mayer expansion 

In this section, we work with a monodisperse system. We note V(u) the interaction 
potential: 

oo if u > 1 

V(u)= , (5) 

otherwise 

The canonical partition function of the system reads: 

Z W = /n rfx k ex p|-E y ( X i- X j- A u)j ( 6 ) 

One would like to study the entropy of the system averaged over disorder: 



S = -\nZ {A} (7) 
In the liquid phase, we can treat this average as an annealed average over the disorder: 



S = S an = - In Z {A} (8) 
So the problem reduces to the study of: 

^=W\I U P ( A im)dA lm f ndxkexp(-£V(xi-xj-A u )) (9) 

Iyyi k \ ij / 

We introduced a 1/JV! prefactor to be able to perform the next step, a Mayer expansion. 
We just have to remember to compensate this factor at the end of the computation. To do 
the Mayer expansion, we have to translate our problem in a grand-canonical formalism. We 
introduce the grand-canonical partition function: 

e = E^w ( 10 ) 

N 
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where z is the activity, related to the chemical potential ji by z = e M , and we rewrite Eq. Q 



as: 



z {A}= f n^n[/(xi-xj)+i] (ii) 

k ij 

where / is an annealed Mayer function: 

7(x - y) = J dAP(A) [exp (-V(x - y - A)) - 1] = - J dAP(A) X (|x - y - A|) (12) 

with x(r) the step function such that x( r ) = 1 if r < 1 and x( r ) = otherwise. We can then 
introduce a diagrammatic representation to express the grand-canonical potential G = In G 
as usuaP^ 

G = In 6 = { connected diagrams } = • + • — •+ j\+ • • • (13) 

where the diagram vertices represent a factor z while the bonds represent the Mayer function 

7- 

B. Mean-field equation of state 

As already noted by KraichnarP, only the first virial correction survives in the mean-field 
limit. The entropy can then be expressed in term of the density field p(x): 

S an [p(x)] = - f dxp(x)[lnp(x)-l] + iy~ dxdyp(x)p(y)7(x-y) + iVlniV (14) 

The IniV contribution to the entropy density reflects the fact that particles are not, in this 
model, truly indistinguishable - just as they are not in a system particles of polydisperse 
sizes. 

The fact that the first terms contribute can be understood by noticing that: 

f(x-y) = -v d /V (15) 



Thus, in the Mayer expansion Eq (13), a diagram with n vertices and m bonds reads 
V n ~ m z n (-Vd) m ■ If m > n, the diagram vanishes in the thermodynamic limit, because this 
requires that several random shifts add to zero. Only the diagrams having m < n contribute, 
and they are the tree ones (having m = n - 1): 

in e = • + •— • + # A. + J~ J + + • • • ( 16 ) 



The next step is to do the Legendre transform^ z p, which leads to: 

In p(x) = In z + ° ® (17) 



where now the ® vertex denotes a factor p(x), and o is 1. Eq. (17) is the saddle point 



equation which minimizes the grand-canonical potential functional G[p(x)], which can be 



translated back into the canonical ensemble, leading to the entropy Eq. (14). 
In the liquid phase with density p, the entropy is thus: 

^l-lnp-^ + lnA (18) 

and the system has a van der Waals equation of state: 

P = P+\v dP 2 (19) 

It is remarkable that this equation of state (ideal gas with first virial correction) is the one 
of hard spheres when d -> oo (see^^). This is a first indication that the limit A -> oo and 
the limit d -> oo are of the same nature, as we shall see below. 



C. Pair correlation function 

As in high-dimensional liquids or in systems with Kac interactions, the very simple form 
of the Mayer expansion allows to compute exactly the two-point correlation functions in the 
liquid phase. A 'naive' pair correlation function is defined as^: 



p(x-y) = < <J(xi -x)5(x 2 -y) > 
N(N-l) 



f l\P(A lm )dA lm 



P 2 N\Z {A} ~> i m 

J Y\ dx k 5(xi - x)5(x 2 - y) exp |- Vfa - xj - Ay) j 



2 _ S\nZ {A} 
-e(xx -x 2 )- 



pN <5e(xi-x 2 ) 



(20) 



From Eq. (14) we then get: 

3(x - y) = 1 + /(x - y) * 1 - v d P(x - y) (21) 
When A = oo the pair correlation function is simply: 

s(x-y) = l (22) 
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This expresses the absence of structure in the system due to the quenched disorder, which 
totally blurs the hard core repulsion when A -> oo. Of course, this does not mean that there 
are no real pair correlations in the system. Indeed, a more interesting quantity is the pair 
correlation function 'seen from one particle': 



ffs(x-y) =<5(x!-x)(5(x2 + Ai2-y) > 
N(N-l) 



., VI7 f U P ( A im)dA lm 

p 2 N\Z {A} J tm 

r / \ (23) 

/ j] dx k <5(xi - x)5(x 2 + A i2 - y) exp I - 2^ ^(x ; - xj - Ay) I 

h \ ij 1 / 

exp(-^(x-y)) 



pA^ 5e(x-y) 



Inserting Eq. (14), we obtain: 

0s(x-y)=exp(-K(x-y)) (24) 

This pair correlation function is identical to the one obtained for a hard sphere system in 
infinite dimensions. 



D. Relation with a Kac potential in the statics of the liquid phase 

The connection between the two ways of approaching a mean-field limit (sending d -*■ oo or 
introducing random shifts with range A -> oo) is to be compared with a similar connection 
in systems with Kac type potentials^. We can push forward this analogy by an explicit 
mapping of our model onto a system with a Kac potential, valid (only) for static quantities 
in the low density liquid phase, when A is large but finite. 



In this case, the Fourier transform of the Mayer function defined by Eq. (12) is the 
product of a Bessel function having a range ~ 1 with the Fourier transform P of the shifts 
distribution, having a range ~ 1/A: 

JAk) 

f(k) = -(2nr"P x (k)^ [fr (25) 
Thus, as long as A » 1, we have: 

f(r) ^ -v d P x (r) = -r d K{r x r) (26) 
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with K a short ranged bounded positive function. This function obviously has a range 7 ~ A, 



as long as A » 1. Thus, (14) reads: 



S m [p{x)] = - f dxp(x)[lnp(x)-l]-^- J dxdyp(x)p(y)K( 7 - 1 |x-y|) + iVlniV (27) 

This is exactly the functional obtained in the so-called mean-field approximation for a 
potential K (with temperature (3 = 1), and it is exact in the Kac limit 7 -*■ ocP^. Indeed, 



the naive pair correlation function Eq. (21) is equivalent to the one obtained for a Kac 
potential. 

This mapping concerns the statics of the liquid phase, and does not mean that the dynam- 
ics of our model has anything to do with the one of a Kac model, even in the liquid phase. 
As a matter of fact, in the following sections we present results showing that a dynamical 
glass transition occurs in our model with a Gaussian distribution for the shifts, whereas the 
related Kac model, which is the Gaussian Core ModeP^ becomes in the high density limit 
an ideal ga d 32 ' 33 1 . 



E. Avoiding the Kirkwood instability 



If A is finite but large, we can take Eq. (14) as a good approximation of the liquid 



phase, but we have to be careful in our choice of the random shift distribution P(A). It 
is well known that the mean-field entropy functional for Kac models can be unstable above 
a given density (the so-called Kirkwood instability^) towards a phase with spatial density 
modulations. This can be seen via a linear stability analysis of Eq. ( Jl4[ ). If we perturb the 
uniform liquid phase solution p(x) = p with a small oscillatory term p e (x) = (l+eexp(ikx))p, 
we get: 

e 2 Np [ 5 2 S 1 

s«[p B (x)] = s a M + — tf [ Mx) ;; (y) j (k) 

with: 



s 2 s n 



5p(x)My) 



~5(x-y) + /(x-y) 



(28) 
(29) 



From Eq. (|28|), it is clear that the uniform liquid phase is stable only if: 

ft 2 9 



TF 



00 



- + /»<o 
p 



Vk 



(30) 



5p(x)5p(y) 

If /(k) takes positive values for some wave vector k, there is a value of p above which the 



condition Eq. (30) is not fulfilled: this is the Kirkwood instability. The only way to avoid 
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FIG. 3. Modulated phase due to a Kirkwood instability. Here the random shift distribution is flat 
within a sphere with radius A = 2. 



this transition is to have a Mayer function with a negative Fourier transform: 

/(k) < Vk (31) 

Since: 

J 4k) 

f(k) = -(2n)^P x (k)^ r (32) 

we have to tune the distribution of the random shifts to get the desired property. If A » 1, 
the range of Pa (~ A -1 ) is much smaller than the range of Jd, therefore -/(k) = VdP\(k) and 



taking P\(k) > is enough to ensure condition (31). In this work, we will take a Gaussian 
distribution for the shifts: 

p ^ A) -(ix)-M-w) (33) 

As an example of the Kirkwood instability, we show in Fig. [3] a dense configuration of the 
random shift model with a flat density of shifts. 



F. Corrections beyond Mean-field: the role of high dimensionality. 

We can consider corrections to the mean-field equation of state. As we shall see, for finite 
A they vanish with dimensionality as oc \~ d . In the Mayer expansion of the entropy, the 
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dominant correction will come from diagrams with m = n, which are the ring diagrams: 

S an>x = NhiN- J dxp(x)[lnp(x)-l]+®-® + g 2 !8 + ^ + ^% + ... (34) 



The resummation of these diagrams has been done by Montroll and Mayei® 

2 r n ../_, \w>) v ' 



San,X 1 , l ll 
= 1 - Inp OVd + -- — — 



and gives, for A » 1 



KPa(A;)] ; 



pv d P x (k) 



+ In 



p 2 / 31 (p) + lniV 



where we have introduced the A-independent factor: 

[v d Pi(k)] a 



lab 



(p)-J 



VLA^&k- 



[l+pVdP^k)} 



The equation of state is then: 



1 2 P 3 

P = P+ 2 P Vd+ (2^W 



(35) 



(36) 



(37) 



(38) 



This expression gives a quantitative estimation of how the approximation becomes better 
at higher dimensions, through the factor X d . Already in d = 3, the ring corrections are less 
than 1% at <fi = 1 for a range as small f as A = 1. When d > 3, one needs to take a very small 
range of random shifts to feel any finite dimensional effect. 



G. Estimation of the density as P -> oo 

The model with A -> oo has a maximal packing density that diverges in the thermo- 
dynamic limit, a rather awkward property from the thermodynamic point of view. When 
A -> oo we can derive an estimation of the maximum density of the random shift model. To 
do this, we consider the following algorithm: 

i) Starting a configuration with spheres in a volume V, generate A^ random shifts, 

ii) Locate positions where one can add, without any overlap, a new sphere interacting with 
other spheres via the A^ random shifts, 

in) If such a position exists, add this (N + l)th sphere, and go back to step 1. 

13 



1 ' — 
10 



100 1000 
N 



10000 



FIG. 4. Maximum density 4> ma x obtained by the algorithm described in the main text, showing 



the lniV dependence of Eq. (41 ) 



Up to which density such an algorithm can work? When A -> oo, the (N + l)th sphere 
will see the other spheres as having completely random positions. Then the probability for a 
position x to satisfy the hard sphere constraints imposed by the N other spheres is simply: 



- W (x) = (l-^) ~-e-^-e-* (39) 



Then the available volume to place the (N + l)th sphere will be Ve~^. The algorithm 
will not be able to find a solution once the available volume is of the order of Vd- If we note 
4>max the volume fraction for which this happens, we have: 

e -^m« „ V A (40) 



which gives: 



In — = lniV-m(lniV) + m(m(lniV))... (41) 



We implemented this algorithm for the random shift model in d = 1. Results are presented 
in Fig. |4j showing the In iV dependence of 4> ma x- 
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III. DYNAMICS IN THE LIQUID PHASE 

A. A warming up exercise: computation of the canonical partition function 

The essential steps followed to compute the exact dynamic equation are similar to the 
ones performed in making an equilibrium computation in the canonical ensemble. It is thus 



instructive to understand them in this context first. Starting again from Eq. (11), and 
introducing the density p(x) by: 

p(x) = £5(x- Xi ) (42) 



with Xi the coordinate of the i-th particle, we can rewrite -Z{a} a s: 

Z^ = fdx l f @[p(x)]8Lx)-ZKx-Xi)\exp[lf dxdyp(x)p(y)ln[l + 7(xi-x j )] 

(43) 

Now we can exponentiate the 5 constraint using a second field p, and integrate over the Xi's: 

Z^ } = J ^[p(x)Mp(x)] 

expjz J dxp(x)p(x) + iVln J dxe~ i/3 ( x ) + - J dxdyp(x)p(y) In [l + /(x ; - Xj)] 

2 (44) 

Again noticing that: 

7(x-y) = -^ (45) 

we can expand the logarithm to get: 

Z^ } = f ^[p(x)Mp(x)] 

r r r 1 r 1 (46) 

expjz J dxp(x)p(x)+iVln J dxe"^ + - J dxdyp(x)p(y)/(x- y)| 

This integral may be evaluated by saddle point with respect to the fields p and p (note that 
each term in the exponential is of order N, including the last one, due to the value of /). 
This gives: 

p-ip(x) 

p(x) = N- 



J dxe-v( x ) (47) 
p(x) =i J dyp(y)J(x-y) 

Thus, the logarithm of the partition function can be written in terms of p as: 

lnZ^ = - f dxp(x)lnp(x) + ^ dxdyp(x)p(y)7(x - y) + NlnN (48) 
Below, we shall follow the same steps, but with trajectories playing the role of particles. 
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B. Derivation of the exact dynamics as a partition function of trajectories 



We wish to study the dynamics of the model in the mean-field limit, starting from the 
Langevin equation: 



Mt) = - E ViV(xi(t) - Xj(t) - Aij) + r]i(t) + hi(t), 



(49) 



where we have introduced an external field hj(i), which acts as a source term, and will later 
be set to zero. The vector rji(t) is a Gaussian white noise with variance 2T, where T is the 
temperature, that we will keep explicit for all the derivation of the dynamics equations: 



< T7i(t)ffe-(f ) >= 2Td Sij 5(t-t') 



(50) 



We denote P({x }, {x T }, r) the probability of having particle i in x° at time t = and 
in xj at time £ = r, in the absence of external field h. We may express this as a sum over 
paths using the Martin-Siggia-Rose/ DeDominicis- Jensen^ formalism. Averaging over the 
noise, this gives: 



P({x°},K},r) = / rjdAyP(A^) 



/{x(0)}={xO} 



with 



{*(*)}] = Tdt 
Jo 



E(x i (t)x l (i) + x l (t)x,(i)) 

- £ (xi(t) Vi^Cxi - Xj - A,,)) + x j (£)V i V(x, - x< - A.,-*)) 

z <ij> 



= £ $ [Xi (t ) , X< (t)] + £ [X, (t), Xi (t) , x, (t ) , Xj (t )] 
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(52) 



where the last two lines define $ and Wij. Integrations in Eq. (51 ) of the variables Xj(t) are 
along the imaginary axis. The quantity P[{xj(t)}, {xj(i)}] = exp [-5[{xj(£)}, (xj(£)}]] may 
seem mysterious, because of the variables Xj(t), which do not have an immediate physical 
meaning. In order to understand them, we consider the probability of paths in the presence 
of external fields (hj(i)}: 



P[{x,(t)}, {x t (t)}, {h,(t)}] = exp {-<S[{x,(t)}, {*,(*)}] - x,(t)h,(t)} , (53) 
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Integrating over the 'hat' variables, we find that: 

p[{x,(f)},{h,(t)}] = f nu[Ht)]P[{xi(t)},{Ht)},{Mt)}] 

r 1 ( 54 ) 

= / u^t)V dt)hi(t) n{^(t)},{ut)}] 

J i 

In other words, the 'hat' variables {xj(t)} are the Fourier-transform variables of the fields: 
the probability of a trajectory ({xj(t)}, (xj(i)}) in the absence of field, is the Fourier trans- 
form of the corresponding physical probability (xj(t)} in the presence of an external field 

{Mt)}- 

The functional formalism casts the dynamical problem into a form that resembles a 
partition function, but with one-dimensional objects (the trajectories) replacing the point 
particles. We may exploit the analogy to repeat the canonical computation in the preceding 
subsection. In particular, we may define the 'dynamical Mayer function' as: 

l + / d [x i (*),x i (0,x J -(t),x j (*)] = y~dA u P(A, J )exp{-^[x l (t),x i (t),x j (t),%(t)]} (55) 

and introduce the dynamical 'partition function': 

Z^ T = f ndx°dx[P {A} ({x },{x^},T) 

r 1 - (56) 

= / n( D N(0,x i (0]e- $[xi(t) ^ (t)] )n( 1 + /i^(0,x i (t),x,(0,x j W]) 

-J A A A 



In the above expression the path integral is now performed over paths with free boundary 
conditions. The initial conditions are weighted with a flat distribution. With these small 



reinterpretations, Z N<T looks like the (static) partition function of iV 'polymers' (x(t),x(t)), 
in an external field interacting via a potential In [l + / d ]. 
Next we introduce a density field: 

p[x(f), x(f )] = Y, 8[*(t) - Xi(t)]6[*(t) - Ut)] (57) 

i 

Note that 5 is here a Dirac function in the sense of trajectories, i.e. a product of ordinary 
deltas, one for each time. Inserting this field in the partition function we obtain: 

n e -»[-i(*)A(*)]=exp[- f D[x(i) ) i(0Mx(t),i(0]*W0,i(0]] 1 (58) 
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and: 



11 ( 1 + /i x * (* ) » x < (*) > x i (* ) > x i (* )] ) = ex P 



^ D[x(t),x(t)]D[y(t),y(t)] 



p[x(t),x(t)]p[y(t),y(t)]ln(l + / (i [x(t),x(t),y(t),y(*)]) 



(59) 



Note that the right-hand side contains a self-interaction term that was not present on the 
left-hand side, but this term is negligible compared to the interparticle interactions. 



As in the static case, looking more closely to the function f d defined by Eq. (55), we see 



that in the integral over disorder, the integrand is 1 if the two trajectories do not interact. 
These trajectories do interact if the random shift is able to bring them close to one another. 
If the trajectories explore a finite volume during between times and r (ie r is not too 
large, at least much smaller than the ergodic time), this is only possible for a finite volume 
of integration on the random shift. As the distribution of the shifts is -P(A) = 1/V, this 
means that the function f d is of order T/V, where T is the typical volume covered by a 
trajectory during an time interval r. Then, just as in the static calculation, we may use 



that ln(l + f d ) ~ f d . Now, imposing the condition Eq. (57) via 
5 



p[x(t), x(t)] - Y, 6 - MtMHt) - Ht)] 



f D[p[x(t),x(t)]] 



3 i/D[x((),x(()>[x(i)i(()]p[x(i),x(i)]-v[x(t),x(i)]5:,i[x(t)-x,(()]J[x(i)-x,(()] 



(60) 



and integrating over Xj(i)'s and Xj(t)'s, we get, for a system at equilibrium at time (we 
dropped the time dependence of the paths to simplify the notation): 



Z NtT = J D[p]D[p]exp^i J D[x,x]p[x,x]p[x,x] + Ann J D[x, x ]e^ [x '* ] 

-J D[x,x]p[x,x]$[x,x] + D[x,x]D[y,y]p[x,x]p[y,y]7 d [x,x,y,x] 



(61) 



For the mean-field dynamics we can take the saddle-point with respect to p and p of the last 
equation. This reads: 

e -ip[x,x] 



P[X '" ] =jV /D[x,x]e-^' 
p[x,x] =-z$[x,x] + z J D[y,y]p[y,y]7 d [x,x,y,y]. 

which gives a closed equation on the density of paths p[x,x]: 

p[x, x] = ^ e _ *[ x '^ + ^ D [ y '^ p [ y '^'' ;d [ x ' :S; ' y '^ 



(62) 



(63) 
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where J\f ensures that the density is normalized to N. 

Reinserting this equation in the partition function, we get the functional S = - In Z, which 
reads: 

S[p] = f D[x,x]p[x,x]ln[p[x,x]] + f D[x, x]p[x, x]$[x, x] + S int [p] - Nln N 

J 1 r J ~ (64) 

s m t[p\ = - g J D [ x > x ] D [y>y]p[ x > x ]p[y>y]/d[ x > x ,y>y] 



Eq. (63) is the counterpart of the saddle point equation (47) obtained in the previous subsec- 
tion for the average density. As happens often in this kind of problems, even the mean-field 
equations are hard to solve, in this case because the complexity of an object like p[x, x] 
makes the problem intractable. But, just as in equilibrium calculations, we can try, as a 



further approximation, to find extrema of the free energy Eq. (63) in a well chosen restricted 
subspace of all possible p[x, x]. The next subsection is devoted to the search of such an 
ansatz. 



C. An approximation in terms of two-point functions 

In terms of physical quantities, the simplest trial form for the probabilities is to propose 
a Gaussian form: 

lnP[{x(t)},{h(t)}] = 

-- J dt a dt b {A a (t a , t b )h(t a )h(t b ) + A 2 (t a , t h )h(t a )x(t b ) + A 3 (t a , t 6 )x(t a )x(t 6 )} (65) 

This leads, equivalently, to proposing an ansatz that is Gaussian in the {x,x} variables. This 
ansatz has to be invariant with respect to a time-independent translation of the trajectory, 
which may be imposed, in terms of the x, x variables, by the form: 

P[{x(t)},{x(t)}] = y~dx exp dt a dt b {B(t a ,t b )±(t a )±(t b ) 

+R(t b ,t a )£(t b )(x(t a ) -x) + R(t a ,t b )Z(t a )(x(t b ) -x) + D(t a ,t b )(x(t a ) -x)(x(t 6 ) -x)} 

(66) 

where the integration over x is over the whole volume and implements the translational 
invariance of the ansatz, i.e. that the quadratic form has a zero-mode in the translations. 
(A similar strategy has been previously used in a static replica ansatz, and it is the same idea 
as the Hill- Wheeler integral that imposes rotational invariance in nuclear theory). It will 
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turn out that, in equilibrium, the ansatz satisfies the causality and fluctuation-dissipation 
relations : 

D(t,t') =0 

1 8 ~ 67 

R(t,t') = ---B(t-t)e(t-t) 

The calculation is cumbersome, and we leave it for Appendix B. The result is expressed 
in terms of the two-time correlation function: 

5(i-0 = ^I<(x a (t)-x Q (t')) 2 > (68) 

where the average is over the Langevin noise and the initial conditions. It may be written 
in the (superficially) Mode-Coupling-like format: 



dB(t) 

dt a 



J dt'E R [B](t-t')B(t') + 2T (69) 



with t > 0. This form is quite general, what defines our equation is the kernel which in 
Mode-Coupling equations is a simple function of B, and here is computed as follows. Define 
first the response function R(t), which satisfies the Fluctuation-dissipation theorem: 

fl(f) = -i^e(-t) (70) 

Then the kernel is: 

Z R (t a -t b ) = | J At a ,dt b , R-\t a -t a ,) ^{*(t a ,)*(t b ,)) int IT\t b , -t b ) (71) 

where R' 1 is the inverse through convolution of R: 

J dt' R(t-t')R- 1 (t') = 5(t) (72) 

and may be obtained from R easily with Laplace transforms. Self-consistency is imposed by 
the definition 

R = Rr 1 

r (73) 
B(t) = J dt" df R- 1 (t-t")B(t n 

The trajectories that contribute to the averages (•)mt are those in which two particles enter 
at any intermediate time within the interaction range, as depicted in figure [5j otherwise their 
contribution vanishes, as in a static Mayer expansion. 
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FIG. 5. The trajectories contributing to a Mayer diagram are those that come into interaction 
range at some time. 

IV. NUMERICAL RESULTS FROM A = oo TO A = 

In order to investigate the system beyond mean-field, we studied the finite A case by nu- 
merical simulations. We used two different Monte-Carlo algorithms, an isobaric-isothermal 
one for the study of the equation of state, and an isovolumic-isothermal one for the study 
of the dynamics close to the glass transition. We worked with systems of 864 particles. The 
system we looked at for the bidisperse case is a 50:50 mixture of particles of diameter 1 and 
1.4, a common choice for a three dimensional bidisperse glass former®^. 

We equilibrate the system at increasing pressures by annealing simulations. We check 
carefully that equilibrium is reached at every studied pressure by performing annealings 
with several compression rates, and ensuring that they give the same values for a given 
observable, e. g. density. Close to the glass transition, we also check that no aging is visible 
in the system. 

A. Simulations in the mean-field case 

In the A -> oo limit, equilibrating the system can be much easier, thanks to the 'planting' 
technique 1 38 1 3 ^ for mean-field problems with quenched random variables, if the annealed free 
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energy is exact (InZ = InZ), one can create an instance in thermal equilibrium by taking 
a random set of variables (here particle positions), and looking for a disorder (here random 
shifts) compatible with this set, ie leading to if ({x}, {A}) = 0. One can show that, within 
mean-field, this creates no bias in the measure. In our model, this means that we can 



generate safely an equilibrium instance as long as the annealed entropy Eq. (14) is valid, 
ie as long as the liquid is the equilibrium phase. Finding an equilibrium instance in this 
ensemble requires at most a few tens of seconds for a system with 1000 particles on a desktop 
computer. This of course, is very useful as we can get easily equilibrated configurations with 
densities close and even above the dynamic glass transition density, as long as we do not 
reach the Kauzmann transition point, if it exists. 



B. Equation of state 

Already at the level of the equation of state, we can notice that the system behaves in a 
mean- field way even for random shifts as small as A = 1. As shown in Fig. [6j all the distance 
between the mean-field and the 3d hard sphere equation of state is covered by systems with 
< A < 1. The inclusion of ring diagrams does not show any noticeable difference with the 
leading order for A > 1, but it takes the agreement down to ~ A = 0.5, which is far beyond 
its expected domain of applicability (A » 1). 

The monodisperse hard-sphere system without disorder undergoes a first order transition 
towards a crystal, as does the system with small values of A. We have checked that nothing 
occurs in the equation of state for systems with shifts as small as A = 0.15. For smaller 
shifts, a weak first order transition is visible, which increases when we decrease A. 



C. Pair correlation function 



We compare the pair correlation function obtained by MC simulations with analytical 



results of section [TTJ In Fig. [7] we show that for <7s(x - y) (Eq. (23)), the analytical form 
derived for large shifts is verified. For finite A, the usual structure of the 3d hard-sphere 
model appears gradually as we reduce A. 
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FIG. 6. The equation of state of the model for different values of the shift range A. Solid lines are the 
analytical equation of state for the mean-field case (for which it is exact), and for the monodisperse 
3d hard-sphere system (the Carnahan-Starling approximation). Points are simulation results. The 
equation of state sticks to the mean- field value for shifts larger than A = 1, except at high pressures 
where a glass transition prevents equilibration. 
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FIG. 7. Pair correlation function <?s(x-y) (Eq. (23)) of a system of N - 864 particles with several 
values of A. 
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V. DYNAMIC GLASS TRANSITION 



In this section, we show that in d = 3, we observe numerically the presence of a dynamic 
glass transition at finite pressure, and we study some dynamical properties close to the 
transition. A numerical simulation cannot exclude that the transition pressure scales with 
the system size (for example as P c ~ In AT), but we have given arguments in the previous 
sections that favor the finite- P c scenario. Note that, on the contrary, we expect that the 
Kauzmann pressure Pk might well scale as logiV, but we shall not try to prove this in this 
paper. 

Our model allows to study various interesting aspects of the dynamic glass transition. 
First, we can follow the qualitative evolution of the glass transition from mean- field to finite- 
dimensions. In particular we ask if we can find in the mean-field limit generic features of the 
mode-coupling theory, and try to see where these features break down as we move away from 
mean-field. Secondly, as we can generate equilibrium mean-field configurations at pressures 
larger than the dynamic glass transition, we can numerically access any desired property in 
the region between dynamic and equilibrium transitions P c < P < Pk (but only at infinite 
A). 

A. Dynamic glass transition 

As is well known, mean field approximations give a dynamic transition (e.g. the mode- 
coupling transition) which is in fact avoided, thanks to activated processes. Here the situa- 
tion is conceptually more clear: for A = oo we expect activated processes to be absent, and for 
finite A we expect them to destroy any pure dynamic transition. For sufficiently large A, one 
may expect that the trace of the dynamic transition is quite clear. The (avoided) dynamic 
transition pressure will still have a dependence on the range, since even the non-activated 
dynamics depends on the A. 

To specify this glass transition we look at three standard quantities: relaxation time, 
diffusion coefficient and dynamic susceptibility Xi(t)- First we look at the relaxation time 
of the system obtained from a two-time correlation function: 

1 N d 

C (9> t )=^ZZ C0S («(*?(«) - x? (*))) (74) 

<J iv i=lo=l 

The curves shown is this article were obtained with q = ir. The relaxation time is defined by 
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FIG. 8. Left Two-time correlation Eq. (74) for the mean- field model for densities between = 0.531 
and (p = 1-7478. There is a plateau appearing upon compression, and the time needed to escape 
this plateau increases strongly with pressure, signaling the existence of a dynamic glass transition. 
Right: Same data rescaled by the relaxation time r a . For high densities, time-density superposition 
holds. 

the time rescaling that gives the best time-pressure superposition in the a regime as shown 
in Fig. U 

The relaxation time shows a super- Arrhenius dependence with pressure (Fig. [9]), and 
seems to diverge for a finite pressure (or density) value, which is characteristic of a fragile 
glass former. This is independent of the range of the disorder. Varying A from oo to 0.2 
increases the glass transition pressure by at most 50%, while the transition density <pd 
drops by almost a factor 2. This is to be related to the large shift in the equation of state 
as we go from a mean-field to a 3d system, as shown in the previous section . 



The same divergence is found for the diffusion coefficient D. In Fig. [TOj we plot the 
relation between D and r a , showing a weak violation of the Stokes-Einstein relation, which is 
commonly observed for supercooled liquids. Note that this curve shows no dependence on A, 
and extrapolates smoothly to the mean-field limit, which is quite remarkable, as the violation 
of Stokes-Einstein relation is often explained by the existence of dynamic heterogeneities, 
which depend on A. 

A natural quantity to study in a fragile glass former is the Xi(t) susceptibility, a measure 
of dynamic heterogeneities. It is defined as the variance of the correlation C: 



X 4(t) = N(<C(q,ty>-<C(q,t) > 2 ) 
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FIG. 9. Left: Relaxation time function of the pressure for several values of A. The 

relaxation is clearly super- Arrhenius, and there is a divergence at a finite pressure for A = oo. For 
finite A, the (apparent) dynamic transition pressure (the analogue of the Mode-coupling transition) 
slowly decreases when A decreases. Right: Relaxation time r a as a function of the density. The 
glass transition density shows a much broader change with the range of the disorder. 
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FIG. 10. Left Diffusion coefficient as a function of the relaxation time for several values of A. 
The diffusion coefficient indicates a glass transition at the same density than the relaxation time. 
Right Evolution of the product Dr a with density. The Stokes-Einstein relation predicts a constant 



For fragile glass former, this quantity exhibits a peak on a time scale ~ r a with a height 
that sharply increases when approaching the glass transition from the liquid (low pressure) 
side. For pressures above the dynamic glass transition (if there is one), we expect x*(t) 
to saturate to a plateau for a time ~ rg, the characteristic time of the /3- relaxation, with a 
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FIG. 11. Left: Dynamic susceptibility Xi(t) (75) for the mean-field model at several densities. 



The peak position and value grow rapidly on approaching to the glass transition. Right: x<i(i) 
with different values of the range of random shifts A, for equal relaxation times. Going from A = oo 
to A = 1 does not affect the dynamic heterogeneities, while for smaller values of A (closer to finite 
a dimensional system), heterogeneities increase strongly. 

height that increases closer to the dynamic transition. When A is finite, we cannot reach the 
glass region, as it is impossible to equilibrate a system with an infinite relaxation time r a . 
However, when A -»• oo, we can use the 'planted-configuration' ensemble mentioned above 
to generate an equilibrium configuration even at a pressure larger than the dynamic glass 
transition pressure. 



A typical result in the mean-field limit of our model for x&{t) is shown in Fig. 11 The 
two expected behaviors - growth from above and from below - are observed, indicating that 
we crossed a dynamic glass transition. The location of this qualitative change coincides with 
the point where the relaxation time seems to diverge, and also to the analytic estimate of 
the transition pressure we give below. 

The dependence of x*(t) 011 ^ ( & t a given relaxation time) again shows little difference 



between the A = 1 and the mean-field cases (see Fig. 11). Further from mean-field, the 



height of Xi(t) depends on A, indicating a clear enhancement of the heterogeneities when 
we get close to the 3d system. This is consistent with the idea that a long ranged disorder 
correlates large regions and thus inhibits heterogeneities occurring on smaller scales. 

We can infer the location of the transition by considering the divergence of the peak value 
Xl of the dynamic susceptibility. This leads to the results compatible with those obtained 
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FIG. 12. Peak of the dynamic susceptibility xl as a function of the density for several values of A. 
This quantity shows a slight difference between A = 1 and the mean-field model. 



from the relaxation time. We plotted the density dependence of the peak value xl 111 Fig- 12 



Note that in the mean- field case, the divergence can be observed on both sides of the glass 
transition. 

It is striking that, for all these quantities, the mean-field behavior seems to be quite close 
to the 3d hard-sphere model. There is no dramatic change between A -> oo and A = 1. This 
is also probably true for other choices of the potential V, suggesting that one can create 
very simply a mean-field caricature of any finite-dimensional glass former by adding random 
shifts with a range of the order of the range of the potential. 



B. The approach to the dynamic transition 

The mode-coupling approximation predicts that the timescale r a diverges algebraically 
with the distance to the glass transition density <p c with an exponent -7. 




In real life, the mode-coupling transition becomes at best a crossover, and the question is 
to what extent should one believe, and in what temperature-pressure range, extrapolations 
within mode-coupling functional forms. In our case, we also expect a divergence in the limit 
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FIG. 13. Left: Peak of the dynamic susceptibility xl as a function of the distance to the dynamic 
glass transition for different values of A. The MCT-like scaling (— — —\ is clearly visible. Right: 
xl as a function of r a in the same regime. The power law relation obtained in MCT also seems to 
hold here for every values of A. 



A = oo. What is interesting about this model, is that we may make A gradually smaller and 
follow the transition as it becomes a crossover, and keep track on the interpolations as they 
become less and less obvious, right down to the original particle model. 



In order to test the behavior (76), we have to fit our simulation data with two free 
parameters <p c and 7. In practice, this is a delicate task as one needs to have a relaxation 
time running over many decades to be able to chose unambiguously the couple {0 C , 7} (see for 
example®'). A way to help this procedure is to look at the four-point dynamic susceptibility 

((f) (j) \ _ 
J . This gives 

us an independent measure of the 7 exponent, using the relation \\ ~ r <^ 7 ^ there is a region 
where mode-coupling scaling holds. We can then look for the pair of parameters {4> c ,l} 
giving the best fit for both measures. 



The results are shown in Fig 13 and Fig. 14 We get an excellent agreement with power 



law divergence on 6 decades of relaxation time for the mean-field model, as expected. For 
small values of A, however, we observe a deviation to the power law scaling when we get close 
to the transition, presumably due to activated processes, a fact that is well attested for a 3d 
binary hard sphere system^. This seems to confirm the idea that mean field theories give 
correct qualitative features for the relaxation time, but only over a limited range of density 
not too close to the glass transition. Quite interestingly, intermediate values of A show a 
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FIG. 14. r a as a function of for <j) c giving the broadest range of densities verifying MCT-like 

behaviour. The agreement is perfect for values of A down to 0.5, whereas for A = 0.2, the dynamic 
transition seems to be missed. 



good power-law scaling up to A = 1, and for A > 1 we were not able to observe any evidence 
of activation within the range of relaxation times that are reachable with our simulations. 

One would like to take these simulations all the way down to A = 0. Unfortunately, this 
is not possible with a monodisperse system, but one can use a bidisperse system and follow 



the same steps. This is what is done in Fig. 15 Here we see that the effect is more clear: 



in term of relaxation time, the algebraic divergence is followed on two decades when A = 0, 



but extends to more than three decades when A = 1. In figure 16 we show the behavior of 
the extrapolated dynamical transition density and the 7 exponent as a function of A. We 
observe that as A is lowered, these values continuously approach the known values for 3d 
systems, both in the monodiperse^ and in the bidisperse case^. 

The relaxation curves for the correlations contain information beyond that of the 
timescale r a . The shape of the relaxation curves does depend on the range of the disorder, 



as is presented on Fig. 17 with the /3-relaxation part becoming slower when A decreases, 
both in the case of monodisperse and of bidisperse systems. There is no visible difference 
between A = 1 and A = 00, yet another indication that for this range of disorder, the system 
behaves in a mean-field manner. For small A, we do not observe a complete separation 
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FIG. 15. Left: r a as a function of for <p c giving the broadest range of densities verifying 

9c 

MCT-like behaviour. The agreement is perfect for values of A down to 0.5, whereas for A = 0.2, the 
dynamic transition seems to be missed. 
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FIG. 16. The values of the parameters of the MCT-like divergence, for different values of A, in 
the monodisperse (left) and bidisperse (right) cases. These quantities show a continuous behavior 
from the 3d, A = system to the mean field A = oo one. 

between a and /3-relaxations, and the plateau is not well defined. These features however 
should strongly depend on the microscopic dynamics, and it is known that Monte-Carlo 
dynamics leads to a longer (3- relaxation than molecular dynamics (see for example^). 

Remarkably, the a-relaxation shape is almost insensitive to the value of A in the monodis- 
perse case. We can just barely notice an increase of the slope when A gets smaller, a feature 
that appears much more clearly in the binary mixture case, as we can access lower A values. 
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FIG. 17. Relaxation curves C{q,t) with different values of the range of random shifts A, for equal 
relaxation times. Left: monodisperse case. The main evolution while tuning A is a lengthening 
the /3-relaxation. Right: Bidisperse case. The evolution is emphasized here. The shape of the 
a-relaxation depends varies strongly between the mean-field case and small values of A. 

C. Small or large exponent? 



We wish to stress again that our goal is this work is to perform a data analysis of the 
same type as performed with MCT, and in particular we wanted to have an exponent 7 that 
fits the divergence of r a on the broadest possible range and satisfies the relation between xl 
and r a . It is clear that by relaxing the constraint given by xl > we can fit t ne divergence of 
T a on a broader range, and possibly on the whole available range, just by taking an exponent 
of the order of four, and a larger density for the divergence point (see for instance®'). 

Indeed, if one looks at the relation xl versus r a very close to the transition (something 
we have not done here but can be found in Fig. 3 of Brambilla et alW^), one can fit the 
relation with a larger 7 on a restricted range. 

Then, one can legitimately think that the exponent 7 ~ 4.5 we find in the mean-field limit 
is also the one of the 3c? hard sphere system, as long as one accepts to relax the constraint 
on the (mode-coupling inspired) relation of xl versus r a . Therefore, what our work proves 
is that, provided we restrict ourselves to orthodox MCT fitting, the glass transition is more 
and more mode-coupling-like when we approach mean-field and the part of mode-coupling 
divergence which is observed in finite dimension is somehow a shadow of the mean-field one. 
What is appealing in this point of view is that the analysis gives for the 3d bidisperse case 
an exponent 7 ~ 2.5, which is compatible with what has already been found in experiments 
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FIG. 18. Left: The function defined in Eq.(95), in terms of the "cage size" parameter K, for 



m = 0.96. The solution for the transition pressure is obtained by searching for the point with the 
largest gradient, which is 1/pd- One also obtains the estimated cage size at this point. Right: 



Density profiles of a single particle in its cage, defined by Eq. (77), obtained from simulations (in 



black), and from analytical approach with a Gaussian ansatz (in red). 



and simulations of comparable systems ^ 7 * 44 * 45 ! , and is reasonably close to the real MCT (not 
only its phenomenology) for a monodisperse system®'. 



D. Analytic estimation of the dynamic transition 

One quick way to determine dynamic (mode coupling- like) transition points, is to make a 
static calculation of the equilibrium state, considered as an ensemble of metastable ergodic 
components; or, equivalently, to look for the lowest pressure at which the effective potential 
(free-energy at fixed distance between configurations) still has two minima. 

In Appendix A we do this, leading to the determination in the figure below: The estimated 
transition density is (pa = 2~ d p d = 1.65 (see left-hand side of figure 18) to be compared with 
the value estimated numerically <pd = 1.82. Note that the analytic computation is not exact 
because it relies on a Gaussian approximation. 

The estimated cage size is also consistent with the simulated values. In the right-hand 



side figure 18, we plot the quantity: 



ffoOO = < [( x *~ < >) - r] > 



(77) 



which is as the density profile of a single particle in its cage, and is the long time limit of 
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the so-called van Hove self-correlation function. 



VI. ONSET PRESSURE 

A decade ago, Sastry et al. introduced a temperature scale T on > T g in supercooled 
liquids, the so-called onset temperature^. They noticed that in a Lennard- Jones binary 
mixture, the energy of the inherent structures (the configurations reached after a quench 
at T = 0) associated with equilibrium configurations at a temperature T shows a crossover 
from a roughly constant value above T on to a regime where it decreases when T decreases. 
They argued that this temperature is the one at which the dynamics becomes landscape- 
influenced, with a super-Arrhenius dependence of the relaxation time. It was a few years 
later arguecP^ that there is a connection between T on and the computed mode-coupling 
density temperature T c , opening the possibility that in mean-field systems T on and T g (= T c ) 
might coincide. The question is legitimate, since one expects that for the spherical p-spin 
glass, the two temperatures do indeed coincide. 

In this section, we show that this connection does not exist for our model. Because we 
are dealing with hard spheres, an inherent structure is the (infinite pressure) configuration 



reached after a rapid compression process, like the one introduced in™ In figure 19 we plot 
the inherent structure density reached after such a process, in terms of the initial density, 
for the model with A = oo . The crossover, corresponding to the 'onset density' (which is the 
relevant quantity for hard spheres, instead of temperature, see previous section) is clearly 
visible. The dynamic transition in this mean-field limit (d = 3) occurs at density <p g = 1.82, 
while we find here <p on ^ 1.15. Because 9 is a well defined quantity here, we see beyond 
doubt that both densities do not coincide. As expected, the crossover density coincides with 



the point where the relaxation (74) starts to show a shoulder. Thus, also in the mean-field 
limit, (p on marks the onset of a qualitative change in the dynamics (see Figj8]), which becomes 
landscape influenced. 



VII. CONCLUSIONS 



We have introduced an approximation scheme for particle systems that is close in spirit 
to the Mode-Coupling approximation, but has the advantage that one may construct a con- 
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FIG. 19. Inherent structure density 4>j as a function of the equilibrium liquid density <p eq . The 
crossover density is clearly different from the dynamic glass transition point, indicated by the 
vertical blue line. 

tinuous range of models, with at one end the original one, and at the other end one for which 
the approximation is exact. The approximation becomes better at higher dimensionality of 
space, everything else remaining equal. 

As it stands, the present scheme is derived from the microscopic model. This was origi- 
nally the case also with the Mode-Coupling equations, although the standard practice has 
become to modify freely the interactions in such a way as to obtain the observed static struc- 
ture factor and transition temperature or pressure. In our case, the analogous procedure 



would be to substitute the true potential by one based on the pair correlation function (23) 



V eff = -T\n[g(r)] (78) 

We have not tried this strategy. 

In this paper we have not discussed the possible static Kauzmann transition, and we 
suspect that it might happen at divergent pressures, of the order of In A, i.e. divergent in 
the mean- field limit. 

By following the model from the original problem to its mean-field limit, we have used 
the present construction to give new arguments on the existence of a vestige of the genuine 
mean-field dynamic transition, following a mode-coupling-like behavior. We also argued 
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conclusively that the 'onset' temperature (or pressure) should not be identified with the 
dynamic transition. More generally, one may follow this strategy to decide whether features 
found in true system that one 'explains' within random first order theory, really extrapolate 
to the corresponding feature in the limit in which the theory is exact. 
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APPENDIX A: AN ESTIMATION OF THE DYNAMICAL TRANSITION 
PRESSURE FROM REPLICAS 

Physical discussion 

The dynamic glass transition is related to the existence of an exponential number of 
amorphous metastable states. Above a given dynamic glass transition density <f>d, the liquid 
phase can be seen as the sum of all these states. 

To derive the properties of the model at high density in order to test the above scenario, 
we will study the partition function of m copies of the original system: 

/m I m \ 

n n dx k° ex p -ee v (x* - v - A u) (79) 
a=l k \ a=l ij ' 

The idea is the following^^ if we force the m copies to be close one another by adding 
a small coupling term between them, we can expect that if we study the system at a density 
higher than 0^, and switch off the coupling after taking the thermodynamic limit, the m 
copies will be confined in the same metastable state. 

This is done in practice by looking at the entropy of the replicated system: 

S = WZ~ = \im-(W l -l) (80) 

where we use the replica trick to average over disorder the logarithm of Z m . Looking more 
closely to Z™, we see that it can be interpreted as the partition function of N 'molecules' 
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Xi = {xj, ...,x™ m } made of nm spheres: 

/r~ ran I ran ^ 

]ldA u P(A u ) / nridx k a exp -EE^-x^-Ay) 
ij J a=\ k \ a=l ij / 

/r~ I van \ 

ridx^n / dA y P(A y )exp -EE^( x i a - x J a - A ij) 
k ij \ a=l ij / 



(81) 



Then, using standard liquid theory techniques, we are able to write the entropy Eq. (80) as 
a functional of the density of molecules p(x). The metastable states of the system are then 
the maxima of this entropy with respect to p(x) when m = 1, which we expect to couple 
replicas. 



Replicated canonical formalism 



We have to study the following partition function Eq. (81 ). Introducing x ; = {xi 1 , Xi" m } 



we can see Eq. (81) as the partition function of a iV molecules interacting via a potential V 
such that: 

// ran \ 
P( A)dAexp -EE V (* a - y a ~ A) (82) 
\ a=l ij ) 

Each 'molecule' Xi is made of n independent sets of m coupled (in the same state) original 
particles. Then, Eq. (81) reads: 

Z»= f ndxin[l + 7 r (xi-Xj)] (83) 

i ij 

We wish to do a Mayer expansion like in the non-replicated case, but we cannot add 



a convenient combinatoric prefactor in (83), as we do not know what it should be a pri- 



ori. However the specific properties of the random-shift model allows us to do a canonical 
treatment of the partition function as we show in the following. Introducing the density of 
molecules p(x) by: 

p(x) = E^( x -Xi) (84) 

i 

we can rewrite as: 

M=f^if D[p(x)]5|p(x)-E5(x-x ? :))exp[^ dxdyp(x)p(y) In [l + / r (x - y)] 

(85) 
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We may now exponentiate the 5 constraint at the cost of adding a second field p, and 
integrating over the Xi's: 

Z% = f D[p(x)]D[p(x)] 

expji J dxp(x)p(x) + Nhx J dxe" 1 ^ + - J dxdyp(x)p(y) In [l + / r (x - y)]| 

(86) 

The next important step is to look more closely at the function / r (x-y), as we did in the 



non-replicated case. In the integral of Eq. (82), the exponential of the potential V vanishes 



1 whenever the two molecules in x and y + A overlap, and gives 1 otherwise. Then we have: 

-nm^<7 ? .(x-y)<"^ (87) 
Thus, by expanding the logarithm, we get: 

Z^ = f D[p(x)]D[p(x)] 

( r r .. 1 r — \ (88) 

expji J dxp(x)p(x) + iVln J dxe^« + - J dxdyp(x)p(y)/ r (x - y) J 

We wish then to evaluate this integral by saddle point with respect to the fields p and p 
(because each term in the exponential is of order N, including the last one, due to the value 
of /). This gives: 

p-ip(x) 

p(x) = N- 



/dxe-v(^) ( 89 ) 
P(x) =i J dyp(y)/ r (x-y) 
Thus, the logarithm of the partition function is: 

\nZ^=-J dxp(x) lnp(x) + i J dxdyp(x)p(y)/(x - y) + N\nN (90) 



We wish to stress the similarity of Eq. (90) with Eq. (14). In fact the average over disorder 
does exactly the same job for the replicated liquid than for the bare non-replicated liquid: 
it disallows 'three- molecule effective interactions', just like in a high-dimensional system. 



Location of the dynamic transition, within the Gaussian ansatz 



We may obtain an approximation for p by assuming it has a Gaussian form. Of course, 



this ansatz will not be a solution of the full saddle point equation Eq. (89), but we can still 
extremize the entropy Eq. (90) within this ansatz. A natural choice is the ansatz^l 

1 



P(x) 



N " 



It, r- 

7=1^ 



dX- 



cxp 



7m 



E 

a=7(m— 1)+1 



(x°-x 7 y 

2K 



(91) 
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Note the similarity with the dynamic ansatz. 



This is a 1-step replica symmetry breaking (1-RSB) ansatz. Injecting Eq. (91) in the 



partition function Eq. (90) leads to integrals exactly similar to the ones one has compute 
in the hard sphere system (without random shifts). This has been done irP^. Following the 
computations along the lines of 5 — one finds: 

S[p(x)] 



(92) 



(93) 



IniV + 1 - lnp - -(1 - m)\n(2nK) + -lnm - -(1 - m) - -I(m, K) 
A 2 2 2 2 

where I(m,K) is the integral: 

I(m, K) = fdX [/ dxdy — exp (^^-) exp (^) X (x - y) 

Then, the saddle point equation on K reads: 

d(m-l) dI(m,K) 
p dlnK 

Thus, there exist metastable states in the liquid phase (which is recovered in the limit 
m -»■ 1) if there is a cage size K which verifies: 



(94) 



1 



lim 



dI(m,K) 



(95) 



p m^id(m-l) <91n.fr 

The dynamic glass transition density is pd, beynd which no solution to this equation can be 
found. It is possible to compute I(m,K) numerically, and we find <pd = ^ d Pd - 1-65 for d = 3 



(see Fig. 18), not far from the value (1.82) found in numerical simulations. The value of K 



at the transition is also comparable with what is found in the simulations (see Fig. 18). 



APPENDIX B : GAUSSIAN ANSATZ FOR THE MEAN-FIELD DYNAMICS 



A. Notation 



The calculation we are going to follow is quite heavy. It may be made somewhat more 
compact, and one may follow the analogy with the static treatment better, by using the 
supersymmetric notation^ (se c 1 54 * 55 *) . One introduces two extra Grassmann variables 6 et 6. 
Denoting a = (t,9,9), the trajectories x(i) and x(t) may be encoded in a superfield ip{a): 



ip(a) = x(t) + 00x(t) 
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(96) 



Defining the operator D a as 

D " = T m§ + <>wm-M < 97 > 

we have: 

$[x,x] = J dadb5(b-a)D a (^(a)-tP(b)) 2 (98) 

Similarly: 

l+7ix,x,x',x'] = J dAP(A)exp|-^ T dtx(t)V x y(x-y-A)) + y(t)V y y(y-x-A)} 
= J dAP(A) exp l~ J da V (^(o) - ip'(a) - A) j 



(99) 



(100) 



P[^] = ~~~r fdxexp[- f dadbB-\a,b) (ip(a) -x)(^(6) -x)l (101) 



The action can be written in the following compact way: 

S[p[V]]= / D^]ln[p[^]] + / D#[^] ^ da^(a)L>^(a) 

By analogy with the statics, one may make a 'Gaussian' ansatz j^. 

V2_ 
V det5(B) 

with 

B(a,o)=0 V a (102) 
In components, the super-correlators read: 

Z3(a, 6) = -B(t a , tfc) - 9 a 9 a R(t b , t a ) - 9 b 9 b R(t a , t b ) + D(t a , t b )9 a 9 a 9 b 9 b 
B-\a, b) = B(t a , t b ) - 9 a 9R(t b , t a ) - 9 b 9 b R(t a , t b ) + D(t a , t b )9 a 9 a 9 b 9 b (103) 

and they are related through: 

J d9 b d9 b dt b B(a, b)B~\b, c) = 5(9 a - 9 c )5(9 a - 6 c )5(t a - t c ) (104) 

This inversion formula may be developed, to obtain the "tilde" variables in terms of the 
ones without tilde. 

In equilibrium, the fluctuation-dissipation theorem implies that B and B~ x take the form: 

B(a, b) = B(t a - t b ) - e a e a R(t b - t a ) - e b e b R{t a - t b ) 

B-\a, b) = B(t a - t b ) - 9 a 9 a R(t b - t a ) - 9 b 9 b R(t a - t b ) (105) 
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where 



R{t) -~T^t ' ~ ~T~dt 



(106) 



where R and R are zero for negative time-differences, due to causality. 
It is easy to show that: 

$[x,x] = - J dx J dad6<y(6-a)D a (V»(a)-x)(V»(6)-x) 
The super-correlator 23 is given by: 

B(a,b) = J D^(a)-V>(&)) 2 p[V>] 



(107) 



(108) 



Inserting this ansatz in Eq. (100), we get: 



5(B) = --TrlnB + J dadb S(a - b)D a B(a,b) - S int (B) 



(109) 



with 



^(B)=^/D[^]/dxdy 

exp[- J dad6B-Ha,6)[(^(a)-x)(^(6)-x) + (^>)-y)(^(6)-y)]] fM^'] 

(110) 

The integrand is invariant with respect to independent translations of x and y, because they 
may be absorbed into the shift A, so we may write: 

Sint(B) = de*Bf D[ ^' ] 6XP [ _ / dad&S_1 (a,6) + 7MV>'] 

(111) 

The saddle point equation gives: 



= 



5S 



5B(a,b) 2 



1 B-\a,b) + 5(a -b)D a 



6B(a,b) 



with: 



SSint 1 a _i 



-B-\a, b) - [B- 1 ® (ip(a')iP(b')) int ® B' 1 ] (a, b) 



and: 



5B(a,6) 2 

Making the convolution product with B(6, c), we obtain: 

= L> a B(a,6) + y dc£(a,c)B(c,6), 
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(112) 

(113) 
(114) 

(115) 



where £ is given by 

E(a, b) = -2 [B- 1 ® (^(a>(&')U ® B" 1 ] (a, 6), (116) 



also of the form (103): 



E(a, 6) = £ B (t a - t 6 ) - 6 a 9 a X R (t b - t a ) - d b d b Y, R {t a - t b ) (117) 
T,b and £# satisfy also a fluctuation-dissipation relation: 

Because correlation and response satisfy a fluctuation-dissipation relation, we may write 
everything exclusively in terms of correlations: 

dB(t a - t b ) 



dt a 



-TR(t b -t a ) + J dt c ^ R {t a -t c )B{t c -t b )-^c{t b -t c )B{t b -t c )t x , (119) 



with: 

d 



E«(ta -h) = - J dt^dty R-\t a - t a ,) ^-(x(t a ,)x(40>mt RT^tv - t b ) (120) 
Given that, as t a ->■ £?(t a - t 6 ) ~ 2T\t a - t b \, we have that: 

dB{ lf b) = -TR(t b , t a ) + f dt c £ R (t a - f e )fl(t c - 1 6 ) + 2T 

OTa (121) 

2T =-Y, B {t b -t c )B{t b -t c )\ h oo , 



which is the result (69) 
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